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Analysis of excess wings in broadband dielectric spectroscopy data of glass forming materials is found to 
provide evidence for anomalous time evolutions and fractional semigroups. Solutions of fractional evolution 
equations in frequency space are used to fit dielectric spectroscopy data of glass forming materials with a range 
between 4 and 10 decades in frequency. We show that with only three parameters (two relaxation times plus one 
exponent) excellent fits can be obtained for 5-methyl-2-hexanol and for methyl-m-toluate over up to 7 decades. 
The traditional Havriliak-Negami fit with three parameters (two exponents and one relaxation time) fits only 
4-5 decades. Using a second exponent, as in Havriliak-Negami fits, the a-peak and the excess wing can be 
modeled perfectly with our theory for up to 10 decades for all materials at all temperatures considered here. 
Traditionally this can only be accomplished by combining two Havriliak-Negami functions with 6 parameters. 
The temperature dependent relaxation times are fitted with the Vogel-Tammann-Fulcher relation which pro- 
vides the corresponding Vogel-Fulcher temperatures. The relaxation times turn out to obey almost perfectly the 
Vogel-Tammann-Fulcher law. Finally we report new and computable expressions of time dependent relaxation 
functions corresponding to the frequency dependent dielectric susceptibilities. 
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I. INTRODUCTION 



Many physical properties of glass forming liquids (e.g. their 
viscosity) vary dramatically (often over 15 or more decades) 
within a narrow temperature interval^. This phenomenon is 
the glass transition. The change of physical properties dur- 
ing the glass transition has not yet been fully understood and 
remains a subject of intense investigations 2 -^. 

In this work we study the glass transition by observing the 
dielectric susceptibility. The dielectric susceptibility quanti- 
fies the response of permanent and induced dipoles to an ap- 
plied frequency dependent electric field. The dielectric loss 
(resp. imaginary part of the complex dielectric susceptibil- 
ity) typically shows a temperature dependent maximum, the 
a-peak, at low frequencies. It is followed at higher frequen- 
cies by a so called excess wing 13 . This excess wing has not 
yet been understood nor has it been described by any model 
with less then 4 fit parameters^. Existing theories, such as the 
mode coupling theory (see 1415 and references therein), do not 
allow to fit the excess wing. Traditional phenomenological fits 
of the excess wing employ a superposition of two Havriliak- 
Negami function s 5,13 ' 16 , and they need 7 fit parameters to fit a 
range of 10 decades. 

The aim of this work is to provide fit functions for ex- 
cess wings with only three (four) parameters (one (two) ex- 
ponents) and two relaxation times) obtained from the previ- 
ously introduced method of fractional time evolutio n 17 ' 18 , and 
to apply them to experimental data exhibiting a clear excess 
wing on a frequency range as broad as possible. Our fitting 
functions need only 3 (model A) or 4 (model B) parameters, 
which is a significant improvement compared to 6 parameters 
for the superposition of the Havriliak-Negami and Cole-Cole 
expression presently used. We study the glass forming mate- 
rials 5-methyl-2-hexanol 6 , glycerol 8 and methyl-m-toluate 12 . 



II. CLASSICAL RELAXATION MODELS 

The Debye relaxation model describes the electric relax- 
ation of dipoles after switching an applied electric field 1 —. The 
normalized relaxation function fit), which corresponds to the 
polarization, obeys the Debye law 



(1) 



T- + l)f(t) = 0, 



with the relaxation time r and initial condition /(0) = 1. The 
response function, i.e. the dynamical dielectric susceptibility, 
X is related to the relaxation function viaii 



x(t) 



dt 



fit)- 



(2) 



In the following discussions we focus on Laplace transformed 
quantities. We use the Laplace transformation of f(t) 



JZ{f{t)}{u) = 



e- ut f(t) dt, 



(3) 



where u = iv and v is the frequency. Rewriting equation © 
in frequency space and utilizing /(0) = 1 leads to 



(4) 



Xiu) = JS?{x(t)}(u) 

= l-uJ?{f(t)}(u) 

1 



1 + UT 



the well known Debye susceptibility. 

In experiments one measures not the normalized quantity 
but instead 



e(u) = (so - £oo)x(«) + £c 



(5) 
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where £o an d £00 are the dynamical susceptibilities at low, 
respectively high frequencies. 

The Debye model is not able to describe the experimental 
data well, because experimental relaxation peaks are broader 
and asymmetric. For this reason other fitting functions 
were proposed like the Cole-Cole 20 , Cole-Davidson— and 
Havriliak-Negami 16 expressions, whose normalized forms 
have typically 2 or 3 parameters (see table H). They were in- 
troduced purely phenomenologically to fit the data. This can 
be considered as a drawback. These functions with three pa- 
rameters are able to fit the data over a range of at most 5 
decades (Havriliak-Negami). Several copies are commonly 
superposed to fit a broader range, e.g. Havriliak-Negami plus 
Cole-Cole, which would result in 6 fit parameters. 



III. FRACTIONAL RELAXATION MODELS 

In this work we use a generalized form of the Debye relax- 
ation model in equation (fTJ. It is based on the theory of frac- 
tional time evolutions for macroscopic states of many body 
systems first proposed in equation (5.5) in 22 and subsequently 
elaborated i n 17 ' 23 ^ 3 - 2 -. As discussed in 17 28 composite fractional 
time evolutions are expected near the glass transition. Such 
time evolutions give rise to generalized Debye laws of the 
form of model A: 



( n D + r 2 tt D a + 1) f(t) = 



or model B: 



(nD 



1 D c 



2 D C 



!)/(*) =0, 



(6) 



(7) 



where the parameters obey < a, a%, 0,2 < 1, a% > a 2 and 
the relaxation times n, T2 > are positive. Here the symbols 
tiD + t 2 q D q , respectively nD + r^D" 1 + r 2 Q2 D Q2 are the 
infinitesimal generators of composite fractional semigroups 
with D" being a generalized fractional Riemann-Liouville 
derivative of order a and almost any type 27 - 33 . If D" rep- 
resents a classical fractional Riemann-Liouville derivative of 
order v then its definition reads (with v £ R + ) 

D*7(*) = DMl"/(i) (8) 
t 

= f^y DM /(* -zr-'fio^ (9) 

fi + u=\v], i >0, 



TABLE I: List of traditional fit functions for dielectric spectroscopy 
data of glass forming materials. 
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Cole-Cole 


1/(1 + (utT) 


2 


Cole-Davidson 


1/(1 + ut)" 


2 


Havriliak-Negami 


1/(1 + («t)T 
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where \v\ is the smallest integer greater or equal v, T the 
gamma function and D M = d M /dt M . 

The Laplace transform of the fractional Riemann-Liouville 
derivative is 34 



J2f{D"/(t)}(«) 



u"Jgf{/(t)}(«) 



With these definitions the Laplace transformation of equa- 
tions © and (O gives with relation (01 the normalized dielec- 
tric susceptibilities of model A 



Xa(u) 



1+T 9 C 



and model B 

Xb(u) 



T\U ■ 



1 



(ID 



(12) 



These results apply also for other types of generalized 
Riemann-Liouville fractional derivatives introduced i n 17 ' 28 . 

The functions from equations (fTTT l and (fT2l are used to fit 
the dielectric spectroscopy data of 5-methyl-2-hexanol, glyc- 
erol and methyl-m-toluate. Real and imaginary part are fitted 
simultaneously with the parameters a, a.\, a 2 , T\ and r 2 . 

Additionally we fit the temperature dependent relaxation 
times t\ and t 2 with the Vogel-Tammann-Fulcher function 



r exp 



DT VF 
T-Tsn 



(13) 



where T is the absolute temperature, tq a material parameter, 
D the fragility and TVf the Vogel-Fulcher temperature. The 
fit parameters are to, D and Typ. 



IV. RESULTS 

Model A fits the 5-methyl-2-hexanol (Fig. [TJ and methyl- 
m-toluate (Fig. [2J data remarkably well for 5, respectively 7 
orders of magnitude, where both the a-peak and the excess 
wing can be fitted simultaneously with only three parameters. 
Model A is better suited then the Havriliak-Negami model 
which only fits reasonably well for up to four orders of mag- 
nitudes for these materials. This improvement is due to the 
positive curvature of the function in (fTTt at frequencies above 
the a-peak. Sometimes this curvature poses also the main dif- 
ficulty when fitting with model A. An example is glycerol as 
seen in the upper part of Fig. [3] While it is easy to fit closely 
the the a-peak it is more difficult to simultaneously fit the ex- 
cess wing. 

Model B can fit the data much better then model A, which is 
not a surprise since it comes with one more parameter. Nev- 
ertheless, it is remarkable that it can fit a range of up to 10 
orders of magnitude with little deviation from the data points. 
We believe that this model can be used to fit over some more 
orders of magnitude, but at this time there is no experimental 
data available which covers a broader range. 
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FIG. 1 : (Color online) Simultaneous fits of real and imaginary part 
with model A (upper figure) and model B (lower figure) for 5-methyl- 
2-hexanol at 155.4 K. Both models show an excellent fitting capa- 
bility. The data are from 6 . 



V. TEMPERATURE DEPENDENCE OF THE 
PARAMETERS 

Because the data have been fitted at different temperatures 
we are able to observe the temperature dependence of the fit- 
ting parameters. For n and t 2 we perform Vogel-Tammann- 
Fulcher fits provided by equation ( [T3V Note that in our nota- 
tion e = e' — ie". From the fits we obtain the Vogel-Fulcher 
temperatures TVfi and TyF2 as well as the fragility parame- 
ters D\ and D 2 for the relaxation times t\ and t 2 for model A 
and model B (see Table HI]). 

For all fits we see a temperature dependence of the relax- 
ation times T\ and t 2 (Fig. [4]- Fig. O that follows the Vogel- 
Tammann-Fulcher fitting function remarkably well. The re- 
laxation times also show a clear downward trend as the tem- 
perature increases, which confirms that t, t\ and r 2 are phys- 
ically meaningful and can be interpreted as relaxation times 
even tough they appear with a non-integer power in equations 
£TQ and (O- 

The parameters a, ot\ and a 2 also show a temperature de- 
pendence. In the case of 5-methyl-2-hexanol (Fig.|4|i there is 
an increase of a with temperature until a plateau near a = 1 
is reached. This effect comes from the decreasing slope of the 
excess wing with increasing temperature. In the fitting func- 
tion of model A this behavior can be achieved by increasing 
a. For the same material there is an apparent increase of a 2 
between 154 K (6.49 K" 1 ) and 287 K (3.48 K" 1 ) which has 
the same origin as the increase in a in model A. By increasing 
a 2 the excess wing becomes less steep. The plateau at 190 K 



FIG. 2: (Color online) Simultaneous fits of real and imaginary part 
with model A (upper figure) and model B (lower figure) for methyl- 
m-toluate at 179.2 K. Model A can fit the data over the whole spec- 
tral range of 7 decades, which is more then with Havriliak-Negami 
which uses the same number of fit parameters. With this data we 
obtain the broadest fit with model A. The data are from 12 . 



(5.26 K -1 ) and above comes from the fact that the fits at those 
temperatures are done mainly for the a-peak since the excess 
wing is not visible. 

For glycerol and methyl-m-toluate there is also a clear tem- 
perature dependence of a, a% and a 2 (Fig. [5] and Fig. [6]). 
The trend is however reversed in comparison to 5-methyl-2- 
hexanol. This comes from the increasing slope of the ex- 
cess wing with increasing temperature. This behavior can be 
achieved in the fit functions by decreasing a, respectively a.%. 



VI. REPRESENTATION OF THE SOLUTIONS AS 
FUNCTIONS OF TIME 

In2£ we obtained the analytical solution of a fractional dif- 
ferential equation of rational order, which we use to analyze 
our fitting results for model A and model B. For a general so- 
lution of equations ( fTTT i and ( fl~2b with arbitrary real cti see^. 
The restriction to rational a* is not a drawback, since we can 
approximate a.\ and a 2 by a rational value on a grid between 
and 1 . This number of grid points is chosen to be 20, which 
keeps computation times reasonably limited as the computing 
time increases quadratically with the lowest common denom- 
inator of a with 1 . 

The solution for f(t) for model B is a sum of Mittag-Leffler 
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FIG. 3: (Color online) Simultaneous fits of real and imaginary part 
with model A (upper figure) and model B (lower figure) for glycerol 
at 185 K. While model A is not able to fit the data well, model B still 
gives an excellent fit over 10 decades. The data are from 8 . 
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FIG. 4: (Color online) Temperature dependence of the fitting param- 
eters for 5-methyl-2-hexanol for model A (upper panels) and model 
B (bottom panels) together with Vogel-Tammann-Fulcher fits. 



type functions: 
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cf- fc - 1 i?(-fc/iV,cf;t), (14) 
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FIG. 5: (Color online) Temperature dependence of the fitting param- 
eters for glycerol for model A (upper panels) and model B (bottom 
panels) together with Vogel-Tammann-Fulcher fits. 
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FIG. 6: (Color online) Temperature dependence of the fitting param- 
eters for methyl-m-toluate for model A (upper panels) and model B 
(bottom panels) together with Vogel-Tammann-Fulcher fits. 



where N is the smallest number for which both a% N and CX2N 
are integers. The coefficients Cj are the zeros of the character- 
istic polynomial 



C N +T^C^ N 



T a 2c a 2 N 



1 = 0, 



(15) 
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FIG. 7: (Color online) The solutions of the fractional initial value 
problems $6$ and ((7} with /(0) = 1 using the fit parameters for 
model A and model B for glycerol at 195 K. 



the function E (v, a; t) is defined as 34 

(at) k 



E{v,a\t)=t v Y J 



k=0 



I> + fc + l) 



(16) 



The coefficients Bj are the solutions of the linear system of 
equations 



JV 

E clB k 

k=l 

0<i<N~ aiN - 1 

N 
fc=l 

N - aiN <i < N - a 2 N 



0. 



(17) 



(18) 



N 

E(4 
fc=i 



_ai i-N+ ai N 



_-«2 _*— N+a 2 N 
T 2 c k 



)B k = 0, (19) 



N - a 2 N < i < N - 2. 



This solution is only valid if all the roots Cj of the character- 
istic polynomial in dT3T > are distinct, which is checked in the 
computations. Since the linear system of equations ([T7]>-(fT9l 
is underdetermined we choose one fundamental solution for 
{Bj} and a multiplication factor for fit") such that /(0) = 1. 

The analytical solutions are plotted for glycerol at 195 K 
(Fig-G)- The fitting values for r for model A are n = 4.991 s 



and t 2 = 1.089 s. Both values lie in the time interval where 
the relaxation occurs, which confirms the interpretation of 
these fitting parameters as relaxation times. For model B the 
fitted times are t\ = 9.729 s and t 2 — 0.92 s. So t 2 marks the 
onset of the relaxation and t\ the end. 

We note that the fractional derivatives appearing in the ini- 
tial value problem ((7} can be generalized to fractional deriva- 
tives of arbitrary type j3 introduced in 27 and defined as 

W^f{t) = l( 1 -«(M-")DMl«M-")/(f) i (20) 
< P < 1. 

For the case j3 = 1 it reduces to the Riemann-Liouville 
fractional derivative, while for j3 = to the Caputo-type 
derivative 36 . Because 



D"'PE (fi, a: t) = W-^E (fx, a; t), 
0</3,7<l, H>-1, v>0, 



(21) 



the solution of our initial value problem does not change by 
replacing the Riemann-Liouville fractional derivatives with 
these generalized Riemann-Liouville fractional derivatives of 
type p. 



VII. CONCLUSIONS 

The two fractional relaxation models (model A and model 
B) are shown to fit well dielectric spectroscopy data for vari- 
ous glass forming materials over a range of at least 5 orders of 
magnitude for model A and 10 orders of magnitude for model 
B. For 5-methyl-2-hexanol and methyl-m-toluate model A can 
fit data over a range of 6 and 7 orders of magnitude. This is 
a significant improvement over conventional fitting formulae 
like Havriliak-Negami and Cole-Cole which need to be su- 
perimposed in order to fit data at the same range. Conven- 
tional fitting formulae require 6 parameters in contrast to 3 
for our model A and 4 for our model B. Another advantage 
of the fractional relaxation models is that the fitting formulas 
follow from an underlying general theory based on composi- 
tion of fractional semigroups and their infinitesimal genera- 
tors, which is not the case for Havriliak-Negami or Cole-Cole 
functions. 
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TABLE II: List of the fit parameters Tvf, D and to for various materials. 
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